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Abstract 

Ab  initio  density  functional  theory,  using  the  B3LYP  hybrid  functional  with  all-electron  basis  sets,  has  been  applied 
to  the  adsorption  of  H  on  the  (0001)  surface  of  wurtzite  GaN.  For  bulk  GaN,  good  agreement  is  obtained  with 
photoemission  and  X-ray  emission  data  for  the  valence  band  and  for  the  Ga  3d  and  N  2s  shallow  core  levels.  A  band 
gap  of  Eg  =  4.14  eV  is  computed  vs  the  experimental  value  (at  0  K)  of  3.50  eV.  A  simple  model,  consisting  of  a  (2  x  2) 
structure  with  3/4-monolayer  (ML)  of  adsorbed  H,  is  found  to  yield  a  density  of  states  in  poor  agreement  with  photo¬ 
emission  data  for  H  adsorbed  on  surfaces  prepared  by  ion  bombardment  and  annealing.  A  new  model,  consisting  of 
co-adsorbed  Ga  (1/4  ML)  and  H  (1/2  ML),  is  proposed  to  account  for  these  data. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  chemisorption  of  H  on  GaN  surfaces  is  a 
key  aspect  of  this  technologically  important  mate¬ 
rial.  The  stability  of  GaN  at  high  temperature  is 
strongly  affected  by  the  presence  of  H2  in  the  ambi¬ 
ent  [i]  which,  in  turn,  influences  the  quality  of 
material  grown  by  metalorganic  chemical  vapor 
deposition  (MOCVD).  Also,  H  remaining  from 
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MOCVD  growth  and  trapped  in  the  bulk,  as  well 
as  at  grain  boundaries  and  interfaces,  affects  the 
performance  of  GaN-based  electronic  devices  [2]. 

Adsorbed  H  is  also  of  fundamental  interest  for 
several  reasons.  The  ideal  Ga-polar  (0001)  and  N- 
polar  (OOOl)  surfaces  are  terminated  in  a  lattice 
plane  of  Ga  or  N,  respectively.  On  such  surfaces, 
Ga  (N)  atoms  have  an  electron  population  of  3/4 
|e|  (5/4  |e|)  localized  in  dangling  orbitals  which 
would  make  these  surfaces  metallic.  On  either  sur¬ 
face,  adsorption  of  0H  =  0.75  monolayers  (ML’s) 
of  H  in  a  (2  x  2)  surface  unit  cell  (where  0H  =  1 
means  one  H  per  surface  site)  satisfies  the 
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“electron-counting  rule”  [3].  This  rule  requires,  for 
a  stable  surface,  that  all  bonding  orbitals  and  N 
dangling  orbitals  be  doubly-occupied  and  that  all 
Ga  dangling  orbitals  be  empty.  On  the  Ga-polar 
face  for  example,  the  three  unpaired  Ga  electrons 
per  (2  x  2)  unit  cell  combine  with  the  three  elec¬ 
trons  on  the  H  atoms  to  form  three  Ga-H  bonding 
orbitals.  This  structure,  henceforth  termed 
“(2  x  2)H”,  leaves  the  dangling  orbital  on  the 
fourth  surface  Ga  atom  empty  and  renders  the  sur¬ 
face  semiconducting  as  observed  experimentally 
(see  below).  Other  stabilization  mechanisms  (e.g., 
vacancies  or  adsorbed  Ga  or  N)  are  also  possible 
as  discussed  in  recent  reviews  [4,5].  Several  theoret¬ 
ical  studies  [3-11]  of  H  adsorption  on  polar  GaN 
surfaces  have  been  performed,  and  most  conclude 
that  the  (2  x  2)H  structure  is  the  most  stable  con¬ 
figuration  except  under  growth  conditions  [3,11] 
involving  a  high  pressure  of  H2.  However,  to  our 
knowledge  the  simple  (2  x  2)H  structure  has  not 
yet  been  observed  experimentally.  A  conceptual 
difficulty  with  the  observation  of  this  structure  lies 
in  how  one  would  stabilize  the  bare  surface,  in  the 
experimentally-observed  semiconducting  state,  be¬ 
fore  adsorption  (or  after  desorption)  of  H. 

A  significant  body  of  experimental  data  exists 
for  the  adsorption  of  H  on  GaN  surfaces; 
although,  a  complete  model  has  not  yet  been  pre¬ 
sented.  Several  groups  have  studied  the  vibrational 
spectrum  [12-17],  using  high-resolution  electron 
energy  loss  spectroscopy  (HREELS),  or  the  tem¬ 
perature-programmed  desorption  (TPD)  [18-21] 
of  adsorbed  H.  Others  have  reported  electron  en¬ 
ergy  loss  (ELS)  [22]  or  ultraviolet  photoemission 
spectroscopy  (UPS)  [23-26]  data  showing  H- 
induced  structure,  which  is  difficult  to  interpret 
without  reference  to  a  theoretical  model.  The 
HREELS  and  TPD  studies  find  no  surface  H  prior 
to  chemisorption;  although,  effusion  of  H  from  the 
bulk  of  MOCVD  material  can  be  detected  [18]  at 
high  temperature.  Furthermore,  UPS  data  show 
that  adsorption  of  H  on  the  clean  surface  removes 
surface  states,  which  indicates  that  such  surfaces 
are  not  already  saturated  with  adsorbed  H.  On 
the  other  hand,  ion  scattering  spectroscopy  (ISS) 
studies  [27,28]  of  the  N-polar  surface  detect  about 
0.75  ML  of  H,  even  on  surfaces  not  intentionally 
exposed  to  H  atoms.  The  surface  H  is  thought  to 


arise  by  out-diffusion,  during  high-temperature 
annealing,  of  bulk  H  retained  from  MOCVD 
growth.  It  has  been  shown  [29]  that  the  release  of 
H  from  the  bulk  of  Mg-doped  (p-type)  GaN  is 
rate-limited  by  the  recombinative  desorption  of 
H2  at  the  surface,  rather  than  by  bulk  diffusion 
of  H  atoms,  and  that  the  surface  coverage  is  low 
under  desorption  conditions. 

It  has  been  found  [30-32]  that,  for  low-energy 
(100  eV)  electrons,  the  cross-section  for  electron- 
stimulated  desorption  (ESD)  of  H  from  the  Ga- 
polar  GaN(0001)  surface  is  several  orders  of 
magnitude  larger  than  for  a  Si(l  00)  or  (1 1 1)  sub¬ 
strate.  The  charge  and  kinetic  energy  of  the 
desorbing  H  and  the  dependence  of  the  desorption 
rate  on  incident  electron  energy,  all  of  which  are 
important  in  defining  the  ESD  mechanism  [33], 
are  unknown  at  present,  and  no  corresponding 
data  have  been  reported  for  the  N-polar  surface. 
It  was  also  found  [31]  that  the  reduction  in  the 
ESD  cross-section  upon  substituting  D  for  H  is 
only  a  factor  of  ~3  for  H/GaN(0001)  vs  ~50 
for  H/Si(100).  This  implies  that,  in  the  former 
case,  desorption  occurs  via  a  very  long-lived  ex¬ 
cited  state  and/or  that  the  kinetic  energy  of  the 
desorbing  particle  is  so  high  that  even  the  more 
massive  D  atom  has  sufficient  velocity  to  escape 
before  relaxation  occurs.  This  issue  is  of  practical 
significance  because  of  the  possibility  of  using 
ESD  of  H  for  patterned  etching  of  GaN  [34]. 

In  the  present  work  we  will  attempt  to  under¬ 
stand  the  changes  in  the  electronic  structure  of  the 
Ga-polar  GaN(000 1)  surface  in  response  to  a  satu¬ 
ration  coverage  of  adsorbed  H.  An  essential  first 
step  is  to  define  a  model  for  the  H-covered  surface 
that  is  consistent  with  the  available  data.  Further¬ 
more,  understanding  the  structure  formed  by 
adsorption  of  H  can  provide  insight  into  that  of 
the  clean  surfaces.  The  focus  here  will  be  on  surfaces 
cleaned  by  ion  bombardment  and  annealing  (IBA) 
following  exposure  to  ambient  atmosphere,  for 
which  the  majority  of  the  available  spectroscopic 
data  (UPS,  ELS,  ISS,  etc.)  have  been  obtained. 
On  the  other  hand,  most  structural  studies,  using 
low-energy  electron  diffraction  (LEED)  or  scanning 
tunneling  microscopy  (STM),  have  been  carried  out 
on  samples  grown  by  molecular  beam  epitaxy 
(MBE)  and  maintained  in  situ.  No  additional  sur- 
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face  cleaning,  and  thus  no  ion  bombardment  or 
high-temperature  annealing  in  ultra-high  vacuum 
(UHV),  is  needed  in  this  case.  The  surface  struc¬ 
tures  resulting  from  the  two  preparation  methods 
are  different,  as  will  be  discussed  below.  Most  of 
the  theoretical  work  1  on  polar  surfaces  to  date 
[3-1 1,35-44]  has  addressed  primarily  the  surface  en¬ 
ergy  under  MBE  or  MOCVD  growth  conditions 
and/or  the  stability  of  the  various  structures  ob¬ 
served  in  LEED  or  STM  for  samples  prepared  in 
situ  by  MBE.  Exceptions  are  the  works  of  Wang 
et  al.  [35,36]  and  of  Strasser  et  al.  [37]  which  inter¬ 
pret  angle-resolved  UPS  data. 


2.  Computational  details 

Density  functional  theory  (DFT)  calculations 
were  performed  using  the  CRYSTAL98  suite  of 
programs  [46,47]  and  the  B3LYP  hybrid  func¬ 
tional.  This  approach  has  been  shown  to  give  rea¬ 
sonably  accurate  band  gaps  for  a  wide  range  of 
materials  [48,49]  and  to  yield  reliable  structural  re¬ 
sults  [50,51]  for  the  bulk  and  for  the  clean  and 
adsorbate-covered  polar  surfaces  of  wurtzite 
ZnO,  which  is  isostructural  and  isoelectronic  with 
GaN.  Optimization  of  basis  sets,  bulk-lattice  con¬ 
stants  and  adsorbate  geometries  (see  below)  was 
done  using  the  LoptCG  script  2  which  computes 
gradients  numerically  by  repeated  runs  of  CRYS- 
TAL98  and  then  executes  a  conjugate  gradient 
procedure  to  locate  the  energy  minimum. 

The  all-electron  contracted  Gaussian  basis 
sets  2  for  Ga  and  N  were  derived  from  those  used 
in  previous  [52,53]  studies  of  GaN.  The  Ga  basis 
set  consists  of  6s  +  5p  +  2d  shells  (an  864111/ 
64111/41  set),  and  the  exponents  of  the  outermost 
3sp  and  Id  shells  were  reoptimized  (using  the 
B3LYP  functional)  with  Ga  in  the  bulk  GaN  lat¬ 
tice.  The  results  were  1.8596,  0.7105  and  0.2763 


1  It  should  be  noted  that  in  some  theoretical  work  the  label 
(0001),  and  not  (0001),  is  applied  to  the  N-polar  face,  and 
comparisons  are  made  to  experimental  data  for  which  the 
surface  polarity  was  not  actually  determined.  Caution  is 
required  on  the  part  of  the  reader  to  avoid  confusion. 

2  The  LoptCG  script  and  the  all-electron  Gaussian  basis  sets 
were  obtained  at  http://www.crystal.unito.it/. 


for  the  3sp  shells  and  0.6393  for  the  d  shell.  The 
original  N  basis  set  consists  of  4s  +  3p  shells.  A 
single  d  shell  was  added  (a  7311/311/1  set),  and 
the  exponents  of  the  outermost  2sp  and  d  shells 
were  optimized.  The  values  found  were  0.4339 
and  0.1319  for  the  2sp  shells  and  0.5972  for  the  d 
shell.  The  standard  Pople  6-3 1G**  basis  set  3 
was  used  for  H,  which  consists  of  inner  and  outer 
s-shells  of  three  and  one  Gaussian,  respectively, 
and  a  single  p-shell  for  polarization.  The  use  of 
DFT  in  CRYSTAL98  requires  auxiliary  Gaus¬ 
sian- type  basis  sets  for  fitting  the  exchange-corre¬ 
lation  potential  [54].  For  Ga  and  N,  these  were 
provided  by  R.  Pandey  (personal  communica¬ 
tion).  4  For  H,  the  standard  set  of  12  s-functions, 
with  exponents  from  0.1  to  2000,  was  used. 

Geometry  optimization  for  the  periodic  slab- 
model  surfaces  with  adsorbates  (see  below)  is 
computationally  intensive.  In  order  to  make  the 
problem  tractable  we  have  adopted  the  practice, 
commonly  used  in  molecular  systems,  of  optimiz¬ 
ing  the  structure  at  a  low  level  and  then  determin¬ 
ing  the  properties  via  a  single-point  calculation  at 
a  high  level.  In  the  present  case,  optimization  of 
the  slab-model  surfaces  was  done  at  the  restricted 
Hartree  Fock  (RHF)  level  using  Durand-Barth- 
elat  effective-core  pseudopotentials  (ECP’s)  as 
described  by  Causa  et  al.  [55]  for  bulk  III-V  sem¬ 
iconductors.  These  are  large-core  ECP’s  that  treat 
the  Ga  3d  level  as  part  of  the  core.  The  21G*  basis 
sets  for  Ga  and  N  were  derived  from  those  given 
[55]  for  GaP  and  BN,  respectively,  with  the  expo¬ 
nents  of  the  outer  sp  and  d  shells  reoptimized 


3  The  H-atom  basis  set  was  obtained  from  the  Extensible 
Computational  Chemistry  Environment  Basis  Set  Database, 
Version  11/29/01,  as  developed  and  distributed  by  the  Molec¬ 
ular  Science  Computing  Facility,  Environmental  and  Molecular 
Sciences  Laboratory  which  is  part  of  the  Pacific  Northwest 
Laboratory,  P.O.  Box  999,  Richland,  Washington  99352,  USA. 
Contact  David  Feller  or  Karen  Schuchardt  for  further  infor¬ 
mation.  The  URL  is  http://www.emsl.pnl. gov:2080/forms/ 
basisform.html. 

4  For  Ga:  14  s-type  functions  with  exponents  from  0.1  to 
6000,  2  d-type  with  exponents  of  0.3  and  0.9  and  1  g-type  with 
an  exponent  of  0.3.  For  N:  14  s-type  functions  with  exponents 
from  0.07  to  2000,  6  p-type  with  exponents  from  0.4  to  1.2,  4  d- 
type  with  exponents  from  0.1  to  1.0,  and  1  f-type  and  1  g-type 
with  exponents  of  0.3  and  0.6,  respectively. 
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for  GaN.  The  values  used  were  Ga:  0.191  (sp), 
0.296  (d)  and  N:  0.255  (sp),  0.800  (d).  As  a  test, 
the  bulk-lattice  constants  obtained  at  the  RHF/ 
ECP  level  will  be  compared  with  those  obtained 
at  the  B3LYP  level  using  the  all-electron  basis  sets 
described  above.  In  calculating  energies  for  the 
various  possible  adsorbate  structures,  the  RHF/ 
ECP  energies  were  corrected  for  electron  correla¬ 
tion  a  postiori  [55]  using  the  HF  electron  density 
and  the  Perdew-Wang  generalized  gradient 
approximation  (PW-GGA)  density  functional. 

Because  the  wurtzite  (0001)  and  (OOOT)  sur¬ 
faces  are  polar,  special  procedures  were  needed 
to  assure  self-consistent  field  (SCF)  convergence 
[56]  in  the  slab-model  calculations.  For  each  sur¬ 
face  structure  (see  below)  the  RHF/ECP  geometry 
optimization  began  with  an  SCF  calculation 
using  a  high  degree  of  Fock  mixing  (95%)  [46]. 
The  resulting  Fock  matrix  was  then  used  as  the 
starting  point  for  another  SCF  calculation 
with  no  Fock  mixing,  and  the  converged  output 
from  this  second  SCF  was  then  used  to  initiate 
the  geometry  optimization.  In  the  DFT  slab- 
model  calculations  it  was  not  possible  to  achieve 
SCF  convergence  with  a  Fock  mixing  of  less  than 
95%. 


3.  Results  and  discussion 

3.1.  Bulk  properties 

Since  this  is,  to  our  knowledge,  the  first  DFT 
treatment  of  GaN  using  a  hybrid  functional,  a 


brief  discussion  of  the  results  for  the  bulk  will  be 
given. 

The  wurtzite  lattice  belongs  to  space  group  186 
(P63mc),  characterized  by  lattice  constants  a  and  c 
and  by  u  -  blc  where  b  is  the  Ga-N  bilayer  thick¬ 
ness.  The  bulk  electronic  structure  was  treated  by 
first  optimizing  a ,  c  and  u  at  the  B3LYP  level  with 
all-electron  basis  sets,  and  the  resulting  parameters 
are  compared  with  experiment  [57-59]  in  Table  1. 
These  values,  rather  than  the  experimental  quanti¬ 
ties,  were  used  in  all  subsequent  calculations  for 
reasons  discussed  elsewhere  [60].  In  dealing  with 
surface  structures,  the  RHF/ECP  method  de¬ 
scribed  above  will  be  used  for  geometry  optimiza¬ 
tion;  hence,  Table  1  also  gives  bulk-lattice 
parameters  obtained  by  this  means.  The  RHF/ 
ECP  results  are  seen  to  be  quite  close  to  the  more 
accurate  DFT/all-electron  values. 

Fig.  1  compares  UPS  data  [61],  extending  from 
the  Fermi  level  to  the  Ga  3d  shallow  core  level, 
with  the  computed  density  of  states  (DOS)  for 
the  optimized  lattice.  The  computed  band  gap,  ob¬ 
tained  from  the  energy  separation  between  the  va¬ 
lence  band  maximum  (VBM)  and  the  conduction 
band  minimum  (CBM)  in  the  DOS,  is  Eg  =  4.14 
eV  vs  the  experimental  value  [57]  of  3.50  eV  at  0 
K.  As  noted  elsewhere  [48],  the  B3LYP  functional 
used  here  is  found  to  overestimate  somewhat  the 
band  gaps  of  several  materials.  The  experimental, 
angle-integrated  spectrum  was  recorded  for  an 
atomically-clean  GaN(0001)  surface  in  ultra-high 
vacuum  under  conditions  of  high  surface  sensitiv¬ 
ity  (hv  =  40.8  eV)  and  high  resolution  (~0. 18  eV). 
Hence,  a  band  of  surface  states  [61]  is  seen  near 


Table  1 

Comparison  of  experimental  and  calculated  bulk  wurtzite  GaN  lattice  constants 


Method 

a  (A) 

c(A) 

da 

u 

Expt.a 

3.1885 

5.1850 

1.6262 

Expt.b 

3.190 

5.189 

1.627 

0.377 

Semi-empir.c 

3.1876 

5.1846 

1.6265 

RHF/ECPd 

3.1937 

5.1876 

1.6243 

0.3801 

B3LYP/all-ee 

3.1899 

5.1866 

1.6260 

0.3761 

a  From  Ref.  [57]  for  a  relaxed  p-type  homoepitaxial  film,  at  room  temperature,  with  a  low  free-electron  concentration. 
b  From  Ref.  [58]  for  a  single  crystal. 

c  Estimated  values  for  pure  GaN  based  on  a  theoretical  study  [59]  of  the  effects  of  impurities  on  GaN  lattice  constants. 
d  This  work,  using  RHF  and  large-core  ECP’s  with  21G*  basis  sets  optimized  for  GaN  (see  text). 
e  This  work,  using  DFT  with  the  B3LYP  functional  and  all-electron  (“all-e”)  basis  sets  optimized  for  GaN  (see  text). 
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Fig.  1 .  Comparison  of  (a)  the  calculated  DOS  for  bulk  wurtzite 
GaN  and  (b)  the  experimental  (surface-sensitive)  UPS.  Binding 
energy  is  referenced  to  the  VBM.  Structure  in  the  DOS  in  the 
vicinity  of  the  Ga  3d  is  shown  on  a  compressed  vertical  scale. 
The  “N  2s”  peak,  labelled  as  is  conventional  in  experimental 
studies,  is  discussed  in  the  text. 


the  VBM  that  is  absent  in  the  bulk  DOS.  Neglect¬ 
ing  the  surface  states,  the  agreement  is  quite  good 
for  the  VB.  Similarly  good  agreement  has  been 
found  in  other  calculations  [62-65]  of  the  GaN 
DOS.  The  Mulliken  charge  on  the  Ga  (N)  atom 
is  +(—)  1.955  |e|.  In  the  slab  calculations  described 
below,  these  will  be  taken  as  benchmark  values  for 
determining  when  a  layer  is  sufficiently  far  re¬ 
moved  from  the  surface  to  be  considered  “bulk¬ 
like”. 

The  Ga  3d  width  and  binding  energy  (BE)  differ 
from  experiment.  The  observed  width  is  affected 
by  the  0.45  eV  spin-orbit  splitting  [66]  (typically 
not  resolved  for  GaN)  and  by  residual  disorder 
which  can  increase  line  width  through  inhomoge¬ 
neous  band  bending.  For  a  partially-ionic  semi¬ 
conductor  like  GaN,  zero-point  lattice  vibration 
can  also  significantly  affect  experimental  core-level 
widths  [67].  None  of  these  effects  are  included  in 
the  calculation. 

Fig.  2  shows  the  Ga  3d  and  N  2s  DOS  in  more 
detail.  This  part  of  the  DOS  has  been  extensively 
discussed  elsewhere  [63-65,68,69].  The  near  reso¬ 
nance  between  the  Ga  3d  and  N  2s  orbitals  in 
the  free  atoms  leads  to  covalent  mixing  which  re¬ 
sults  in  bonding  and  anti-bonding  bands,  together 


-18  -16  -14  -12 


Energy  Below  VBM  (eV) 

Fig.  2.  Total  (a)  and  partial  (b,c)  DOS’s  for  the  Ga  3d  and  N  2s 
levels.  The  bonding  character  of  each  band  is  indicated,  and  the 
anti-bonding  band  is  repeated  for  each  trace  at  10-fold 
magnification. 


with  non-bonding  states  which  are  nearly  pure  Ga 
3d.  The  bonding  character  of  the  different  bands  is 
shown  [69]  in  charge  density  plots.  While  the  non¬ 
bonding  and  anti-bonding  bands  are  easily  seen  in 
UPS  (Fig.  1),  the  bonding  band  has  only  recently 
been  observed  [70,71],  in  N  K-shell  X-ray  emission 
spectroscopy  (XES).  The  dipole  selection  rule  re¬ 
quires  that  states  observed  in  N  K-shell  XES  have 
a  finite  N  2p  component.  Orbital  decomposition  of 
the  DOS  (not  shown)  did  indicate  such  a  contribu¬ 
tion  and  that  it  was  weaker  than  the  N  2s 
component. 

The  discrepancy  between  the  DOS  and  the  UPS 
data,  as  regards  the  Ga  3d  and  N  2s  BE’s,  has  been 
explained  by  Fambrecht  et  al.  [64]  in  terms  of  the 
inability  of  DFT,  with  a  local  and  energy-inde¬ 
pendent  exchange-correlation  potential,  to  model 
single-particle  excitations  of  narrow  bands.  How¬ 
ever,  Vogel  et  al.  [72,73]  have  shown  that  the 
use,  within  the  local  density  approximation,  of 
self-interaction-  and  relaxation-corrected  pseudo¬ 
potentials  virtually  eliminates  the  disagreement 
with  experiment  in  the  vicinity  of  the  Ga  3d.  Sim¬ 
ilar  improvements  in  the  BE’s  for  shallow  core 
levels  were  obtained  for  a  wide  range  of  par¬ 
tially-ionic  semiconductors.  However,  the  GaN 
VB  appears  to  be  relatively  unaffected  by  the 
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corrections  [73],  and  the  band  gap  obtained  (4.0 
eV)  was  close  to  that  found  in  the  present  work. 

The  present  B3LYP  results  are  in  somewhat 
better  agreement  with  experiment  than  are  previ¬ 
ous  results  which  were  based  on  fully  ab  initio 
functionals  and  which  also  neglected  the  self-inter- 
action  and  relaxation  corrections.  The  energy 
difference  between  the  non-bonding  and  anti¬ 
bonding  bands,  about  2.7  eV,  is  close  to  the  value 
of  ^3  eV  observed  for  the  “Ga  3d”  and  “N  2s” 
peaks  (Fig.  1).  The  computed  non-bonding  band 
lies  at  15.6  eV  below  the  VBM  vs  the  experimental 
value  [74]  of  17.76  eV,  and  the  bonding  band  falls 
at  1.1  eV  below  the  non-bonding  band,  in  good 
agreement  with  the  experimental  value  [71]  of 
1.2  eV. 

3.2.  Adsorption  of  H 

3.2.1.  The  (2  x  2)  H  surfaces 

The  first  goal  is  to  determine  a  structure  for  the 
H-adsorbed  surfaces  that  will  adequately  model 
the  data  described  above.  Many  theoretical  studies 
(e.g.,  Refs.  [3-11,35-43])  have  shown  that  the 
polar  GaN(000 1)  and  (000  T)  surfaces,  when  com¬ 
pletely  clean  and  free  of  defects,  are  energetically 
unfavorable  in  comparison  to  surfaces  with  either 
vacancies  or  adsorbates.  As  noted  above,  either 
autocompensation  via  vacancy  formation  or  the 
presence  of  adsorbates  is  needed  to  satisfy  the  elec¬ 
tron-counting  rule.  Since  the  main  focus  is  on  H- 
adsorption,  the  bare  ideally-terminated  surfaces 
are  not  analyzed  here. 

We  begin  with  a  brief  examination  of  the  well- 
known  [3,6-11]  (2  x  2)H  surfaces,  described  above, 
to  see  if  they  are  useful  for  the  present  purposes. 
Fig.  3a  shows  the  model  for  the  (0001)-  and 
(000l)-(2  x  2)H  surfaces  based  on  a  two-dimen- 
sionally-periodic  slab  comprised  of  five  Ga-N 
bilayers.  Since  the  present  basis  sets  use  localized 
atomic  orbitals,  rather  than  plane  waves,  periodic¬ 
ity  along  the  surface  normal  is  not  needed,  and  an 
isolated  slab  can  be  used  [50,51].  Previous  calcula¬ 
tions  have  used  slabs  of  four  [3,6,9,35],  five  [8]  or 
eight  [39,41]  bilayers,  and  it  will  be  shown  below 
that  a  five-bilayer  slab  is  sufficiently  thick  to  isolate 
the  two  surfaces  from  each  other.  The  Ga-H  and 
N-H  distances  were  set  at  1.58  and  1.07  A,  respec- 


Fig.  3.  (a)  Slab  model  for  the  (0001)-  and  (0001)-(2x2)H 
surfaces  with  0H  =  0.75  ML.  The  slab  consists  of  5  Ga-N 
bilayers  with  three  of  the  four  atoms  on  either  surface  bonded 
to  H.  (b)  One  of  several  possible  models  for  co-adsorbed  Ga 
and  H  on  the  Ga-polar  face  with  the  Ga  adatom  in  the  T4  site. 
The  N-polar  face  is  shown  H-free  with  a  Ga  adatom  in  the  H3 
site. 

tively,  based  on  the  sums  of  covalent  radii  of  the 
respective  atoms.  A  Ga-H  distance  of  1.59  A 
was  previously  found  [8]  for  the  optimized 
(000 1)-(2  x  2)H  structure.  In  the  present  work 
the  central  (i.e.,  third)  bilayer  was  essentially 
bulk-like  in  that  the  atom  charges  were  ±1.958 
|e|,  virtually  identical  to  the  bulk  value,  noted 
above,  of  1.955  |e|.  In  contrast,  surface  Ga  (N) 
atoms  bonded  to  H  had  charges  of  +1.80  (-1.60) 
|e|.  The  H-atom  charges  were  -0.439  |e|  for  Ga- 
H  and  +0.322  |e|  for  N-H. 

The  DOS  results  (see  below)  for  the  (2  x  2)H 
slab  are  in  qualitative  agreement  with  those  de¬ 
scribed  by  Fritsch  et  al.  [8]  and  may  be  summa¬ 
rized  as  follows.  5  The  Ga-H  bonding  orbitals  on 
the  Ga-polar  face  are  located  in  a  narrow  region 
(~0.5  eV  wide)  near  the  VBM.  On  the  N-polar 
face,  the  N  doubly-occupied  non-bonding  lone- 
pair  orbital  produces  a  narrow  level  in  the  gap,  just 
above  the  VBM,  and  the  N-H  bonding  orbitals 


5  The  band-structure  results  illustrated  in  Ref.  [8]  are  for 
AIN,  but  the  GaN  results  are  similar  and  are  described  in  the 
text. 
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occur  near  the  bottom  of  the  VB  over  a  range  of 
~2.5  eV.  In  the  present  case,  geometry  optimiza¬ 
tion  was  not  done,  but  the  results  concur  with 
those  of  Fritsch  et  al.  [8]  for  which  such  optimiza¬ 
tion  was  performed.  Adsorbate  features  are  very 
weak  in  the  partial  DOS  (PDOS)  for  the  central 
GaN  bilayer.  Together  with  the  bulk-like  charge 
densities  on  the  central  bilayer  (noted  above),  this 
supports  the  conclusion  that  a  five-bilayer  slab  is 
sufficiently  thick  to  isolate  the  two  surfaces  from 
each  other. 

The  VB  DOS’s  obtained  for  the  simple  (000 1)- 
(2x2)H  and  (000l)-(2  x  2)H  surfaces  bear  little 
resemblance  to  UPS  data  reported  to  date  (see  be¬ 
low).  The  H-induced  features  seen  in  UPS  on  the 
Ga-polar  face  [23]  are  broad  and  spread  essentially 
across  the  whole  VB.  On  the  N-polar  face  [26], 
there  is  no  evidence  for  either  N-H  bonding  states 
near  the  VB  minimum  or  an  N  non-bonding  lone- 
pair  orbital  above  the  VBM.  The  latter  state 
occurs  at  ~1.4  eV  below  the  VBM  on  the  clean 
surface  [26,35]  and  is  removed  by  H  adsorption. 
One  is  then  led  to  conclude  that  such  (2  x  2)H 
structures  do  not  adequately  represent  the  experi¬ 
mental  surfaces,  which  have  been  prepared  by 
IB  A  followed  by  exposure  to  H  atoms. 

Another  model  for  H  adsorption  was  consid¬ 
ered;  namely,  one  involving  a  Ga-H-Ga  bridge 
in  which  a  three-center,  two-electron  (3C2E)  bond 
is  formed.  Such  structures  are  familiar  in  the 
chemistry  of  boranes  (B^H^),  and  an  In-H-In 
3C2E  bridge  has  been  observed  [45]  in  the  chemi¬ 
sorption  of  H  on  the  In-polar  InP(OOl)  surface. 
In  that  case,  an  In  atom  on  an  ideally-terminated 
surface  has  two  dangling  orbitals,  each  with  an 
effective  occupancy  of  3/4  |e|.  Electron  counting 
can  then  be  satisfied  by  forming  a  doubly-occupied 
In-In  dimer  bond  with  an  In-H-In  3C2E  bridge 
between  dimers.  However,  an  analogous  structure, 
with  electron-counting  requirements  fully  satisfied, 
cannot  easily  form  on  the  GaN  surfaces  studied 
here. 

Yet  another  model  could  be  formulated  by  add¬ 
ing  only  one  H  atom  to  the  (2  x  2)  unit  cell,  giving 
6>h  -  0.25  ML.  Of  the  resulting  four  “surface” 
electrons  (i.e.,  three  from  the  surface  Ga  atoms 
and  one  from  the  H),  two  would  form  the  Ga-H 
bond  and  two  would  form  a  doubly-occupied 


non-bonding  lone-pair  orbital  on  one  of  the  other 
Ga  sites.  However,  this  structure  violates  the  elec¬ 
tron-counting  rule  described  above  and  has  not 
been  further  investigated  here. 

3.2.2.  Models  for  co-adsorbed  Ga  and  H 

If  adsorbed  H  is  not  directly  and  solely  respon¬ 
sible  for  stabilizing  the  polar  surfaces,  under  the 
conditions  of  interest  here,  then  other  mechanisms 
must  be  considered.  Based  on  previous  theoretical 
work  [3-11,35-43],  Ga  or  N  adatoms  or  vacancies 
are  potential  candidates.  The  field  can  be  nar¬ 
rowed  by  recalling  that  HREELS  data  [12,13]  after 
exposure  of  the  Ga-polar  surface  to  H  atoms  show 
only  Ga-H  vibrations.  This  argues  against  a  signif¬ 
icant  concentration  of  Ga  vacancies  or  N  adatoms 
on  this  surface,  which  would  be  expected  to  lead  to 
the  appearance  of  N-H  modes.  The  temperature  at 
which  air-exposed  GaN  samples  are  typically  an¬ 
nealed  in  UHV,  as  part  of  the  IBA  cleaning  proc¬ 
ess,  is  in  the  800-900  °C  range.  This  is  above  the 
temperature  (780  °C  [75])  for  congruent  decompo¬ 
sition  of  GaN  and  corresponds  to  a  more  rapid 
loss  of  N  than  of  Ga.  Hence,  a  reasonable  surface 
model  in  the  present  case  may  involve  Ga  adatoms 
remaining  from  incomplete  thermal  desorption  of 
the  outermost  Ga-N  bilayer. 

Indirect  evidence  for  the  possible  importance 
of  Ga  adatoms  in  H  adsorption  comes  from 
HREELS  data  [15-17]  for  H  adsorbed  on 
GaN(0001)  which  show  (unlike  the  case  for  the 
Ga-polar  face  [12,13])  both  Ga-H  and  N-H 
stretching  modes.  Unfortunately,  the  HREELS 
resolution  (typically  ~5— 10  meV)  is  not  sufficiently 
high  that  multiple,  distinct  Ga-H  sites  can  be  eas¬ 
ily  distinguished  on  the  Ga-polar  surface.  Further 
evidence  comes  from  the  observation  [26]  in  UPS 
of  a  Ga-derived  surface  state  on  (000T)  surfaces 
prepared  by  IBA.  The  state  is  removed  by  chemi¬ 
sorption  of  H. 

The  adsorption  of  Ga  on  the  Ga-  and  N-polar 
faces  has  been  studied  theoretically  [3-11,35-43], 
mainly  for  conditions  pertaining  to  MBE  growth. 
For  the  (0001)  face  under  Ga-rich  conditions, 
adsorption  of  6>oa  =  0.25  ML  in  a  (2  x  2)  structure 
leads  to  a  stable  surface  (henceforth  termed 
“(2  x  2)Ga”)  with  the  adatom  in  a  T4  site  (directly 
above  an  underlayer  N  atom).  This  structure  is 
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slightly  lower  in  energy  than  one  with  the  adatom 
in  an  H3  site  (directly  above  a  3-fold  hollow  in  the 
N  underlayer).  The  T4-H3  energy  difference  is  pre¬ 
dicted  to  be  small  in  comparison  to  A Hads,  the  heat 
of  adsorption  of  H  (see  below).  The  three  Ga-Ga 
backbonding  orbitals  are  all  doubly-occupied,  and 
the  dangling  orbitals  on  the  fourth  lattice  Ga  and 
on  the  Ga  adatom  are  empty,  which  satisfies  elec¬ 
tron-counting  requirements. 

A  (000 1)-(2  x  2)Ga  structure  has  been  reported 
in  STM  studies  [44,76,77]  following  in  situ  MBE 
growth.  It  has,  however,  been  suggested  [76,78] 
that  this  results  from  a  small  amount  of  co¬ 
adsorbed  As  impurity,  due  to  previous  use  of  the 
MBE  chamber  to  grow  As  compounds.  In  any 
case,  a  well-developed  (2  x  2)  LEED  pattern  is 
not  generally  seen  for  surfaces  prepared  by  IBA 
(Refs.  [29,79]  and  works  cited).  Instead  a  (1  x  1) 
pattern  is  found  (often  accompanied  by  faceting) 
which  suggests  a  poorly-ordered  surface,  possibly 
with  a  random  array  of  Ga  adatoms.  This  is  differ¬ 
ent  from  the  pseudo-(l  x  1)  adlayer  structure, 
involving  two  full  monolayers  of  adsorbed  Ga, 
which  is  seen  under  Ga-rich  MBE  conditions 
[38,39].  For  the  IBA  (1  x  1)  surface,  UPS  data 
show  no  evidence  of  metallic  character  in  the  form 
of  emission  at  the  Fermi  level,  and  a  high  coverage 
of  Ga-Ga  bonds  is  not  detected  in  surface-sensi¬ 
tive  Ga  3d  XPS  (e.g.,  Ref.  [61]).  For  a  monolayer 
of  Ga  on  the  GaN  surface  [80],  the  adlayer  Ga 
3d  appears  as  a  distinct  shoulder  at  ~  1.4  eV  lower 
BE  than  that  of  bulk  GaN,  and  a  Ga  bilayer  re¬ 
sults  in  an  even  larger  chemical  shift.  It  has  been 
shown  [81,82]  that  a  (2x2)  surface  (which  indi¬ 
cates  a  Ga-polar  face  [83])  can  be  obtained  when 
air-exposed  samples  are  cleaned  by  annealing  in 
NH3  vapor  6  (without  ion  bombardment);  how¬ 
ever,  another  study  [84]  found  only  a  (1  x  1)  sur¬ 
face  for  a  similar  treatment. 


6  A  (000 1)-(2  x  2)N  structure,  involving  N  adatoms,  is 
known  to  be  intrinsic  to  the  Ga-polar  face  grown  under  N-rich 
MBE  conditions  (Ref.  [78]  and  works  cited).  It  is  possible  that 
this,  rather  than  Ga  adatoms,  might  be  the  origin  of  the  (2  x  2) 
structures  reported  in  Refs.  [81,82].  However,  for  reasons  given 
above,  it  is  considered  unlikely  that  a  significant  coverage  of  N 
adatoms  occurs  on  the  IBA-cleaned  Ga-polar  surface. 


In  the  following,  we  will  invoke  a  (2  x  2)Ga 
structure  as  a  model  for  the  local  environment  of 
a  Ga  adatom  and  use  this  as  a  starting  point  for 
the  subsequent  adsorption  of  H.  This  is  based  on 
the  theoretical  stability  of  this  structure  and 
on  the  experimental  evidence,  noted  above,  for 
the  importance  of  Ga  adatoms.  We  put  aside  the 
question  of  whether  a  pure  (2  x  2)Ga  phase  exists, 
under  typical  conditions,  as  a  stable,  macroscopic 
entity  on  the  Ga-polar  surface.  The  Ga  adatom 
is  initially  placed  in  a  T4  site;  however,  due  to 
the  small  T4-H3  energy  difference  (see  below),  it 
is  necessary  to  consider  the  possibility  that  H 
adsorption  causes  a  displacement  of  the  adatom 
to  the  H3  site.  It  is  noted  in  this  regard  that  the 
activation  barrier  for  hopping  between  T4  and 
H3  is  calculated  [85]  to  be  0.4  eV. 

Adsorption  of  two  H  atoms  on  a  (2  x  2)Ga  sur¬ 
face  leads  to  structures  like  that  shown  in  Fig.  3b, 
which  continue  to  satisfy  the  electron-counting 
rule.  One  of  the  three  Ga-Ga  adatom  backbonds 
is  broken,  and  there  are  then  six  possible  structures 
as  shown  in  Fig.  4.  Whichever  Ga  is  vacant  (i.e., 
not  bonded  to  H  or  to  the  Ga  adatom)  then  exhib¬ 
its  an  unoccupied  dangling  orbital.  Including  the 
possibility  of  T4  — >  H3  displacement  then  leads  to 
a  total  of  12  inequivalent  structures.  In  the  follow¬ 
ing,  structures  will  be  labelled  “Ga(T4)  +  2H(a)”, 
etc.  to  specify  the  Ga  and  H  sites  with  reference 
to  Fig.  4.  Finally,  in  defining  the  structures  in 
Fig.  4,  no  distinction  has  been  made  between 
inequivalent  Ga  sites  2  and  4  (see  Fig.  4a).  Addi¬ 
tional  unique  structures  could  be  formed  by  mov¬ 
ing  the  adsorbed  H  from  site  4  to  site  2  in  models 
(a)  and  (c).  Since  these  sites  are  related  by  an 
approximate,  local  mirror  plane  (not  a  true  mirror 
plane)  we  expect  these  to  be  energetically  equiva¬ 
lent,  which  will  be  justified  below. 

The  first  step  was  to  obtain  the  geometries  and 
energies  of  the  H-free  (2  x  2)Ga(T4)  and  Ga(H3) 
structures.  These  results  are  needed  in  estimating 
A//acis,  and  they  also  serve  as  a  check  of  the  calcu- 
lational  procedure.  In  either  case,  the  adatom  posi¬ 
tion  was  optimized  while  keeping  the  rest  of  the 
slab  fixed.  The  height  of  the  Ga(T4)  adatom  above 
the  terminating  Ga  plane  was  found  to  be  1.66  A 
vs  previously-reported  values  of  1.78  A  [8],  1.66 
A  [38]  and  1.73  A  [44].  For  Ga(H3)  the  height 
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Fig.  4.  Schematic  diagram  showing  six  possible  structures, 
derived  from  the  (2  x  2)Ga  surface,  after  adsorption  of  two  H 
atoms  per  unit  cell.  The  solid  circles  are  H  atoms,  and  the  solid 
lines  show  backbonds  from  the  Ga  adatom  to  the  surface.  In 
(a),  the  inequivalent  Ga  sites  are  numbered  for  reference  in  the 
text.  The  adatom  is  shown  in  the  T4  site,  and  another  six 
structures  could  be  drawn  with  the  H3  site  occupied.  These  are 
model  structures  only  and  have  not  been  optimized.  Also  shown 
is  a  diagram  for  the  H-free  (2  x  2)Ga  unit  cell  with  the  Ga 
adatom  in  the  T4  site  and  with  the  empty  H3  site  indicated. 


above  the  surface  plane  was  1.72  A  vs  a  previous 
value  [38]  of  1.63  A.  The  H3  site  was  found  to 
be  higher  in  energy  than  the  T4  by  0.69  eV  per 
(2  x  2)  unit  cell  based  on  the  RHF/ECP  energies 
after  correction  for  electron  correlation  as  de¬ 
scribed  above.  This  compares  with  previous  results 
of  0.6  eV  [8]  and  0.12  eV  [38].  A  vanishingly- small 
energy  difference  was  found  in  Ref.  [85].  In  these 
and  in  subsequent  calculations  the  N-polar  face 
was  “passivated”  with  one  Ga  adatom  per  (2  x  2) 
unit  cell  in  the  H3  site  which  has  been  found 
(e.g.,  Ref.  [8])  to  be  the  most  stable  structure  for 
this  face  at  low  Ga  coverage. 

Previous  studies  of  the  (000 1)-(2  x  2)Ga(T4) 
structure  found  that  the  only  relaxation  of  the  sub¬ 
strate  is  a  small  outward  displacement  of  the 


Table  2 

Computed  heats  of  adsorption,  per  (2  x  2)  unit  cell,  for  the 
(2  x  2)Ga  +  2H  co-adsorbed  structures  shown  in  Fig.  4 


Structure*1 

Atfads  (eV) 

Ga(T4)  +  2H(a) 

6.33 

Ga(T4)  +  2H(b) 

6.35 

Ga(T4)  +  2H(c) 

8.30 

Ga(T4)  +  2H(d) 

6.32 

Ga(T4)  +  2H(e) 

6.24 

Ga(T4)  +  2H(f) 

8.43 

Ga(H3)  +  2H(a) 

6.49 

Ga(H3)  +  2H(c)b 

8.05 

Ga(H3)  +  2H(f) 

8.03 

a  The  labelling  of  the  structures  is  described  in  the  text  and 
in  Fig.  4. 

b  For  this  structure  (referring  to  Fig.  4c)  the  H  atom  on  site  1 
is  moved  to  site  3,  and  the  Ga  adatom  backbonds  are  to  sites  1 
and  2. 


outermost  Ga  layer  by  0.13  A  [8]  or  0.10  A  [44], 
relative  to  the  ideally-terminated  position.  An 
even  smaller  relaxation,  0.05  A,  is  found  [8]  for 
the  (2  x  2)H  surface  described  above.  Hence,  in 
the  present  study,  only  the  Ga  and  H  adatom  posi¬ 
tions  were  adjusted  in  the  geometry  optimizations. 

All  six  of  the  Ga(T4)  +  2H  structures  indicated 
in  Fig.  4  were  subjected  to  geometry  optimization 
at  the  RHF/ECP  level.  The  coordinates  of  each  of 
the  three  adsorbed  atoms  (one  Ga  plus  two  H) 
were  optimized  while  atoms  in  the  slab  itself  were 
frozen  in  the  ideal,  bulk-lattice  positions.  For  each 
structure,  A//ads  (Table  2)  was  obtained  using 

E[(2  x  2)Ga(T4)]  +  2  *  E[ H] 

=  E[G a  +  2H]  +  2  *  FZP(Ga-H)  +  A//ads 

where  E[( 2  x  2)Ga(T4)]  and  if[Ga  +  2H]  are, 
respectively,  the  total  RHF/ECP  energies  of  the 
(2  x  2)Ga  surface  before  and  after  H  co-adsorp¬ 
tion.  Both  values  were  corrected  a  postiori  for  elec¬ 
tron  correlation  as  described  above.  E,ZP(Ga-H)  is 
the  zero-point  vibrational  energy  of  an  adsorbed  H 
atom,  and  if[H]  is  the  energy  of  a  free  H  atom. 
Based  on  HREELS  data  [12,13],  EZF(Ga- 
H)  =  0.116  eV.  Any  H-induced  change  in  the 
vibrational  modes  of  the  slab  is  neglected  in  esti¬ 
mating  the  total  EZF.  A  value  of  if[H]  =  -13.558 
eV  was  found  for  the  present  H-atom  basis  set 
which  is  close  to  the  experimental  value  of 
—  13.606  eV.  The  A Hads  computation  uses  two  H 
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atoms,  rather  than  an  H2  molecule,  since  the 
adsorption  experiments  involve  dissociation  of 
H2  by  a  hot  (^2000  K)  tungsten  filament.  Molec¬ 
ular  H2  is  not  known  to  chemisorb  readily  on  the 
GaN(OOOl)  surface  at  room  temperature  under 
UHV  conditions. 

The  results  in  Table  2  indicate  that  Ga(T4)  + 
2H(f)  is  the  most  energetically-favorable  structure. 
The  optimized  interatomic  distances  are  given  in 
Fig.  5a.  Various  aspects  of  the  results  are  notewor¬ 
thy.  Firstly,  adsorption  of  H  on  the  Ga  adatom  is 
unfavorable  in  comparison  to  having  both  H’s  on 
lattice  Ga  sites.  Mulliken  population  analysis  (at 
the  RHF/ECP  level)  shows  that  the  effective 
charge  on  the  Ga  adatom  in  Ga(T4)  +  2H(f)  is 
+0.31  |e|;  whereas,  the  charge  on  the  Ga  atoms 
bonded  to  H  is  +1.07  |e|.  Given  the  H  effective 
charge  of  about  -0.37  |e|  it  is  evident  that  the  lat¬ 
tice  Ga  sites  should  be  favored.  Secondly,  the  vir¬ 
tually  identical  Hads  values  for  (a)  and  (d)  and  for 
(b)  and  (e)  indicate  that,  all  else  being  equal,  it 
makes  little  difference  which  of  the  three  Ga-Ga 
backbonds  is  broken.  This  supports  the  approxi¬ 
mation,  noted  above,  of  treating  Ga  sites  2  and  4 
as  being  equivalent. 

A  third  point  is  that  the  AHa(±s  values  are  quite 
large,  e.g.,  4.2  eV/Ga-H  bond  for  Ga(T4)  +  2H(f) 
vs  4.5  eV  for  the  bond  energy  in  H2  [86].  Strain 


(a)  Ga(T4)+2H(f)  (b)  Ga(T4)+2H(a) 

Fig.  5.  (a)  Schematic  diagram  showing  the  optimized  surface 
structure  for  the  lo  west-energy  model  (see  Table  2), 
Ga(T4)  +  2H(f).  For  clarity,  only  the  uppermost  two  Ga-N 
bilayers  are  shown,  (b)  Similar  to  (a)  but  for  a  structure  with 
one  H  on  the  Ga  adatom,  Ga(T4)  +  2H(a).  The  adatom  Ga-H 
bond  makes  an  angle  of  7.8°  with  the  adatom  Ga-Ga-Ga 
plane. 


in  the  Ga-Ga  backbonds  for  the  H-free 
(2  x  2)Ga  surface  probably  makes  a  significant 
contribution  to  A Hads.  In  structures  with  one  H 
on  the  adatom  the  local  geometry  of  the  Ga  ada¬ 
tom  is  seen  to  be  nearly  planar  (Fig.  5b)  as  ex¬ 
pected  for  3-fold  coordinated  Ga  in  a  free 
molecule  [86].  When  one  of  the  three  backbonds 
is  broken,  the  adatom  is  free  to  relax  away  from 
the  ideal  T4  position  and  into  this  nearly-planar 
configuration. 

The  A Hads  obtained  here  is  not  in  good  agree¬ 
ment  with  the  experimental  value  [29]  of,  at  most, 
2.44  eV/Ga-H  bond  which  is  based  on  the  kinetics 
of  hydrogen  release  at  700-900  °C.  However,  the 
configuration  of  the  H-covered  surface  under  such 
conditions  may  differ  from  that  in  the  0  K  limit 
appropriate  to  the  present  calculations.  For  exam¬ 
ple,  if  desorption  at  ^  700  °C  occurs  not  from 
the  Ga(T4)  +  2H(f)  configuration  but  from  the 
Ga(T4)  +  2H(a)  as  an  intermediate  state,  then  a 
A//acis  of  about  3.12  eV/Ga-H  bond  would  be  ex¬ 
pected  (cf.  Table  2),  in  somewhat  better  agreement 
with  experiment.  Including  the  thermal  excitation 
of  the  Ga-H  vibrational  mode  at  ^  700  °C  would 
further  reduce  the  theoretical  estimate  of  the 
desorption  energy. 

A  final  point  is  that  the  presence  of  adsorbed  H 
appears  to  reduce  the  T4-H3  energy  difference.  As 
noted  above,  the  present  RHF/ECP  result  places 
the  H-free  Ga(H3)  site  0.69  eV  higher  than  the 
Ga(T4);  whereas,  for  example,  the  Ga(H3)  +  2H(f) 
structure  lies  only  0.40  eV  higher  than  the 
Ga(T4)  +  2H(f).  Although  the  effect  of  H  on  the 
activation  barrier  for  T4  H3  hopping  is  un¬ 
known,  the  results  suggest  that  co-adsorption  of 
H  during  MOCVD  might  possibly  promote  two- 
dimensional  growth  by  enhancing  Ga  adatom  sur¬ 
face  diffusion. 

We  have  also  obtained  the  energy  difference, 
AE ,  between  the  Ga(T4)  +  2H(f)  structure  and 
the  (2  x  2)H  surface,  described  above,  after  optim¬ 
izing  the  positions  of  the  three  H  atoms  in  the 
latter  unit  cell.  An  optimized  Ga-H  distance  of 
1.60  A  was  obtained.  The  energy  difference  is  given 
by 

AE  =  E[ Ga(T4)  +  2H(f)]  -  E[( 2  x  2)H] 

-  £(Ga)  +  £(H) 


V.M.  Bermudez  /  Surface  Science  565  (2004)  89-102 


99 


where  if(G a)  and  if(H)  are  the  energies  of  the 
respective  free  atoms  obtained  for  the  basis  sets 
used  in  the  slab  calculations.  The  result  is 
A E  -  —2.5  eV,  i.e.,  the  Ga(T4)  +  2H(f)  is  the  more 
stable  structure. 

3.2. 3.  Electronic  structure  of  co-adsorbed  Ga  and  H 
Having  determined  the  most  energetically-favo- 
rable  structure  we  will  now  examine  its  electronic 
properties  using  DFT  with  the  B3LYP  functional 
and  the  all-electron  basis  sets  described  above. 
Figs.  6  and  7  show  the  valence  band  UPS  differ¬ 
ence  spectrum,  A N(E),  for  adsorbed  H  in  compar- 


-16  -12  -8  -4  0 


Fig.  6.  Comparison  of  UPS  data  for  H/GaN(0001)  with 
calculated  results,  (a)  Hell  UPS  difference  spectrum  (adsorbed 
minus  clean),  Ref.  [30].  The  resolution  is  about  0.8  eV.  (b,c) 
Calculated  H-derived  states  for  the  (2  x  2)H  and 
Ga(T4)  +  2H(f)  surfaces.  All  binding  energies  are  referenced 
to  the  VBM.  For  trace  (a),  the  VBM  was  found  by  linear 
extrapolation  of  the  VB  edge  to  zero  intensity.  For  traces  (b) 
and  (c)  the  VBM  was  determined  from  the  PDOS  of  the  inner¬ 
most  Ga-N  bilayer.  The  vertical  scale  for  (b)  has  been  reduced 
by  a  factor  of  10  relative  to  that  of  (c). 


-16  -12  -8  -4  0 


Binding  Energy  (eV) 

Fig.  7.  Similar  to  Fig.  6  but  showing  a  comparison  of  the  data 
with  the  calculated  PDOS  for  the  four  surface  Ga  atoms  and 
the  one  Ga  adatom  in  the  (2  x  2)  unit  cell  (cf.  Fig.  4).  (a)  Same 
as  in  Fig.  6.  (b,c)  The  Ga(T4)  surface  before  and  after  formation 
of  the  Ga(T4)  +  2H(f)  structure.  Binding  energy  is  referenced  to 
the  VBM.  Si  and  S2  are  surface  states  on  the  bare  (H-free) 
surface  that  are  affected  by  adsorption  of  H. 

ison  with  the  computed  PDOS  derived  from 
adsorbed  H  (Fig.  6)  and  from  surface  and  adatom 
Ga’s  (Fig.  7).  The  A N(E)  data  were  obtained  [30] 
by  subtracting  the  clean-surface  spectrum  from 
that  recorded  after  H  adsorption.  The  clean-sur¬ 
face  spectrum  was  shifted  (to  account  for  a  change 
in  band-bending)  and  scaled  so  as  to  minimize  the 
residual  Ga  3d  intensity  in  A N(E).  As  is  typical  in 
angle-integrated  UPS  studies  of  adsorption  on 
GaN(OOOl)  [23,61,87],  a  “negative”  feature  ap¬ 
pears  in  A N(E)  just  above  the  VBM  due  to  re¬ 
moval  of  surface  states  characteristic  of  the  clean 
surface  prepared  by  IBA. 

Fig.  6b  shows  the  H  PDOS  for  the  (2  x  2)H 
model,  described  above,  which  concurs  with  previ¬ 
ous  results  [6,8].  As  noted  above,  the  H-derived 
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states  on  the  Ga-polar  face  are  concentrated  near 
the  VBM,  in  conflict  with  the  data  which  show 
these  states  to  be  distributed  essentially  through¬ 
out  the  whole  YB.  Much  better  agreement  is  found 
for  the  Ga(T4)  +  2H(f)  structure.  It  is  useful  to 
consider  why  the  H  PDOS  is  so  different  for  the 
two  structures.  For  the  Ga(T4)  +  2H(f),  adsorp¬ 
tion  may  be  visualized  as  occurring  in  three  steps; 
namely,  breaking  a  Ga-Ga  backbond,  redistribu¬ 
tion  of  charge  to  form  two  singly-occupied  Ga 
dangling  orbitals  and,  finally,  formation  of  two 
Ga-H  bonding  orbitals.  The  intermediate  step  is 
a  relaxed  excited  state  of  the  initial  structure,  and 
the  wavefunction  of  this  state  will  be  a  combina¬ 
tion  of  many  occupied  and  virtual  orbitals  of  the 
ground  state.  This  is  reflected  in  the  spread  of 
the  H  contribution  over  a  wide  energy  range  in 
the  final  adsorbed  state.  For  the  (2  x  2)H,  on  the 
other  hand,  the  ground  state  of  the  (hypothetical) 
bare  (2  x  2)  surface  unit  cell  already  exhibits  occu¬ 
pied  Ga  dangling  orbitals,  and  no  extensive  mixing 
of  occupied  and  virtual  orbitals  is  needed  to  “pre¬ 
pare”  the  surface  for  adsorption. 

Fig.  7b  shows  two  Ga-derived  surface  states  on 
the  H-free  surface.  State  Si  lies  below  the  VBM 
and  is  localized  largely  on  the  adatom.  State  S2, 
which  is  derived  from  adatom  Ga-Ga  backbonds, 
lies  in  the  gap  at  ~1.2  eV  higher  energy  than  the 
experimentally-observed  surface  state.  These  re¬ 
sults  are  qualitatively  consistent  with  those  of  Frit- 
sch  et  al.  [8]  which  is,  to  our  knowledge,  the  only 
previous  work  to  present  surface-state  bands  spe¬ 
cifically  for  the  (000 1)-(2  x  2)Ga(T4)  structure.  5 
Their  results  show  an  Si  state  (comprised  of  ada¬ 
tom  s-orbitals)  lying  slightly  above  the  VBM  and 
an  S2  adatom  backbonding  state  lying  somewhat 
farther  into  the  gap.  In  the  present  work,  adsorp¬ 
tion  of  H  (Fig.  7c)  strongly  attenuates  and/or 
broadens  Si.  The  S2  state  shifts  and  broadens 
slightly  in  response  to  H  and  decreases  in  intensity 
when  one  of  the  three  adatom  backbonds  is  bro¬ 
ken.  However,  S2  is  not  completely  eliminated  by 
H  adsorption  which  is  an  apparent  discrepancy 
with  experiment  [25]. 

Comparison  can  also  be  made  to  the  angle- 
resolved  UPS  results  of  Chao  et  al.  [25].  These 
were  obtained  for  a  sample  of  unspecified  polarity 
which,  however,  was  grown  by  MOCVD  on  a  Si- 


polar  6H-SiC(0001)  surface  and  which  we  thus 
assume  to  be  Ga-polar.  Two  non-dispersive,  H- 
sensitive  states  were  observed  at  or  just  below  the 
VBM  on  the  clean  surface  which  could  correspond 
to  Si  and  S2.  A  third,  highly-dispersive  surface 
state  was  also  observed  which  is  not  found  in  the 
present  calculations.  On  the  N-polar  face,  a  similar 
highly-dispersive  surface  state  [26]  has  been  inter¬ 
preted  [35]  in  terms  of  a  (lxl)  Ga  adlayer  with 
interacting  Ga  dangling  orbitals.  Such  a  structure 
is  not  included  in  the  present  model.  Also,  no  H- 
induced  features  were  observed  by  Chao  et  al. 
[25];  however,  A N(E)  spectra  were  not  presented, 
which  could  mean  that  the  very  broad  and 
nearly-structureless  emission  due  to  Ga-H  orbitals 
was  overlooked. 


4.  Summary  and  conclusions 

Ab  initio  computational  methods  have  been  ap¬ 
plied  to  the  adsorption  of  H  on  the  GaN(0001) 
surface.  The  results  are  as  follows. 

(1)  Density  functional  theory  with  the  B3LYP 
hybrid  functional  gives  a  good  description  of 
the  bulk  lattice  and  electronic  structure  of 
wurtzite  GaN.  In  particular,  the  results  for 
the  hybridization  of  the  Ga  3d  and  N  2s  levels 
is  in  somewhat  better  agreement  with  experi¬ 
ment  than  are  previous  results  using  other 
functionals. 

(2)  The  ideal  H-stabilized  (000 1)-(2  x  2)H  struc¬ 
ture,  with  0.75  ML  of  adsorbed  H,  does  not 
account  for  UPS  data  for  H  on  GaN(0001) 
surfaces  cleaned  by  ion  bombardment  fol¬ 
lowed  by  high-temperature  annealing  in 
UHV. 

(3)  A  new  model  is  proposed  which  gives  results 
in  better  agreement  with  UPS  data.  Here, 
0.50  ML  of  H  adsorbs  on  a  (2  x  2)Ga  struc¬ 
ture  with  0.25  ML  of  adsorbed  Ga.  Both  H 
atoms  in  the  (2  x  2)  unit  cell  adsorb  on  sur¬ 
face-plane  Ga  sites,  leaving  the  Ga  adatom 
free  of  H. 

The  model  developed  here  provides  a  reasona¬ 
bly  good  description  of  the  UPS  data  for  IBA  sur- 
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faces  with  adsorbed  H  and  is  consistent  with  other 
results  suggesting  the  presence  of  Ga  adatoms  on 
H-free  surfaces  prepared  by  IBA.  Although  the 
model  provides  a  starting  point  for  understanding 
the  unusually  large  ESD  cross-section  for  H/ 
GaN(000 1),  it  has  not  yet  provided  any  direct  in¬ 
sight  into  the  mechanism  for  this  effect.  For  exam¬ 
ple,  evidence  was  looked  for  (but  not  found)  to 
suggest  a  significantly  higher  degree  of  Ga-H  vs 
Ga-N  bond  co valency.  Such  an  effect,  were  it  to 
exist,  would  enhance  the  rate  for  Auger  decay  of 
a  Ga  3d  core  hole  via  a  process  involving  one  or 
both  electrons  in  the  Ga-H  bond,  relative  to  decay 
through  an  interatomic  process  involving  the  N 
nearest-neighbors  in  the  partially-ionic  GaN 
lattice. 
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